The theory of canonical perturbations appHed to 
attitude dynamics and to the Earth rotation. 

Osculating and nonosculating Andoyer variables. 

Michael Efroimsky 

us Naval Observatory, Washington DC 20392 USA 

^ ! and 

o ; Alberto Escapa 

Departamento de Matematica Aplicada, Universidad de Alicante, Alicante E-03080 Spain 

February 2, 2008 

Abstract 



> 

■ In the method of variation of parameters we express the Cartesian coordinates or the 
T^lj- , Euler angles as functions of the time and six constants. If, under disturbance, we endow 

' the "constants" with time dependence, the perturbed orbital or angular velocity will 

consist of a partial time derivative and a convective term that includes time derivatives 
. of the "constants". The Lagrange constraint, often imposed for convenience, nullifies the 

I convective term and thereby guarantees that the functional dependence of the velocity on 

Qh' the time and "constants" stays unaltered under disturbance. "Constants" satisfying this 

Q ■ constraint are called osculating elements. Otherwise, they are simply termed orbital or 

rotational elements. When the equations for the elements are required to be canonical, it 
^ ' is normally the Delaunay variables that are chosen to be the orbital elements, and it is the 

Andoyer variables that are typically chosen to play the role of rotational elements. (Since 
some of the Andoyer elements are time-dependent even in the unperturbed setting, the 
^ ' role of "constants" is actually played by their initial values.) The Delaunay and Andoyer 

■ sets of variables share a subtle peculiarity: under certain circumstances the standard 
equations render the elements nonosculating. 

In the theory of orbits, the planetary equations yield nonosculating elements when 
perturbations depend on velocities. To keep the elements osculating, the equations must 
be amended with extra terms that are not parts of the disturbing function (Efroimsky 
& Goldreich 2003, 2004; Efroimsky 2005, 2006). It compUcates both the Lagrange- and 
Delaunay-type planetary equations and makes the Delaunay equations noncanonical. 

In attitude dynamics, whenever a perturbation depends upon the angular velocity 
(like a switch to a non-inertial frame), a mere amendment of the Hamiltonian makes the 
equations yield nonosculating Andoyer elements. To make them osculating, extra terms 
should be added to the equations (but then the equations will no longer be canonical). 

Calculations in nonosculating variables are mathematically valid, but their physical 
interpretation is not easy. Nonosculating orbital elements parameterise instantaneous 
conies not tangent to the orbit. (A nonosculating i may differ much from the real incli- 
nation of the orbit, given by the osculating i.) Nonosculating Andoyer elements correctly 
describe perturbed attitude, but their interconnection with the angular velocity is a 
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nontrivial issue. The Kinoshita-Souchay theory tacitly employs nonosculating Andoyer 
elements. For this reason, even though the elements are introduced in a precessing frame, 
they nevertheless return the inertial velocity, not the velocity relative to the precessing 
frame. To amend the Kinoshita-Souchay theory, we derive the precessing-frame-related 
directional angles of the angular velocity relative to the precessing frame. 

The loss of osculation should not necessarily be considered a flaw of the Kinoshita- 
Souchay theory, because in some situations it is the inertial, not the relative, angular 
velocity that is measurable (Schreiber et al. 2004, Petrov 2007). Under these circum- 
stances, Kinoshita's formulae for the angular velocity should be employed (as long as 
they are rightly identified as the formulae for the inertial angular velocity). 

1 The Hamiltonian approach to rotational dynamics 

1.1 Historical preliminaries 

The perturbed rotation of a rigid body has long been among the key topics of both spacecraft 
engineering (Giacaglia & Jefferys 1971; Zanardi & Vilhena de Moraes 1999) and planetary 
astronomy (Kinoshita 1977; Laskar & Robutel 1993; Touma & Wisdom 1994; Mysen 2004, 
2006). While free spin (the Euler-Poinsot problem) permits an analytic solution in terms of 
the elhptic Jacobi functions, perturbed motion typically requires numerical treatment, though 
sometimes it can be dealt with by analytical means (like, for example, in Kinoshita 1977). 
Perturbation may come from a physical torque, or from an inertial torque caused by the frame 
noninertiality, or from nonrigidity (Getino & Ferrandiz 1990; Escapa, Getino & Ferrandiz 2001, 
2002). The free-spin Hamiltonian, expressed through the Euler angles and their conjugate 
momenta, is independent of one of the angles, which reveals an internal symmetry of the 
problem. In fact, this problem possesses an even richer symmetry (Deprit & Elipe 1993), 
whose existence indicates that the unperturbed Euler-Poinsot dynamics can be reduced to 
one degree of freedom. The possibility of such reduction is not readily apparent and can 
be seen only under certain choices of variables. These variables, in analogy with the orbital 
mechanics, are called rotational elements. It is convenient to treat the forced-rotation case as 
a perturbation expressed through those elements. 

The Andoyer variables are often chosen as rotational elements (Andoyer 1923, Giacaglia & 
Jefferys 1971, Kinoshita 1972), though other sets of canonical elements have appeared in the 
literature (Richelot 1850; Serret 1866; Peale 1973, 1976; Deprit & Elipe 1993; Fukushima & 
Ishizaki 1994)0. After a transition to rotational elements is performed within the undisturbed 
Euler-Poinsot setting, the next step is to extend this method to a forced-rotation case. To this 
end, one will have to express the torques via the elements. On completion of the integration, 
one will have to return back from the elements to the original, measurable, quantities - i.e., to 
the Euler angles and their time derivatives. 

1.2 The Kinoshita-Souchay theory of rigid-Earth rotation 

The Hamiltonian approach to spin dynamics has found its most important application in the 
theory of Earth rotation. A cornerstone work on this topic was carried out by Kinoshita (1977) 
who switched from the Euler angles defining the Earth orientation to the Andoyer variables, 

^ Some authors use the term "Serret-Andoyer elements." This is not correct, because the set of elements 
introduced by Richelot (1850) and Serret (1862) differs from the one employed by Andoyer (1923). 
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and treated their dynamics by means of the Hori (1966) and Deprit(1969) techniqueo Then he 
translated the results of this development back into the language of Euler's angles and provided 
the precessional and nutation spectrum. Later his approach was extended to a much higher 
precision by Kinoshita & Souchay (1990) and Souchay, Losley, Kinoshita & Folgueira (1999). 



1.3 Subtle points 

When one is interested only in the orientation of the rotator, it is sufficient to have expressions 
for the Euler angles as functions of the elements. However, when one needs to know also the in- 
stantaneous angular velocity, one needs expressions for the Euler angles' time derivatives. This 
poses the following question: if we write down the expressions for the Euler angles' 
derivatives via the canonical elements in the free-spin case, will these expressions 
stay valid under perturbation? In the parlance of orbital mechanics, this question may be 
formulated like this: are the canonical elements always osculating? As we shall demonstrate 
below, under angular-velocity-dependent disturbances the condition of osculation is incompat- 
ible with that of canonicity, and therefore expression of the angular velocity via the canonical 
elements will, under such types of perturbations, become nontrivial. 

In 2004 the question acquired a special relevance to the Earth-rotation theory. While the 
thitherto available observations referred to the orientation of the Earth figure (Kinoshita et 
al. 1978), a technique based on ring laser gyroscope provided a direct measurement of the 
instantaneous angular velocity of the Earth relative to an inertial frame (Schreiber et al. 2004, 
Petrov 2007). 

Normally, rotational elements are chosen to have evident physical interpretation. For ex- 
ample, the Andoyer variable G coincides with the absolute value of the body's spin angular 
momentum, while two other variables, H and L , are chosen to coincide, correspondingly, with 
the Z-component of the angular momentum in the inertial frame, and with its z-component 
in the body frame. The other Andoyer elements, g , I , h , too, bear some evident meaning. 
Hence another important question: will the canonical rotational elements preserve their 
simple physical meaning also under disturbance? 



2 The canonical perturbation theory 
in orbital and attitude dynamics 

2.1 Kepler and Euler 

In orbital dynamics, a Keplerian conic, emerging as an undisturbed two-body orbit, is regarded 
to be a "simple motion," so that all the other available motions are conveniently considered as 
distortions of such conies, distortions implemented through endowing the orbital constants Cj 
with their own time dependence. Points of the orbit can be contributed by the "simple curves" 
either in a nonosculating fashion, as in Fig. 1, or in the osculating manner, as in Fig. 2. 

The disturbances, causing the evolution of the motion from one instantaneous conic to 
another, are the primary's nonsphericity, the gravitational pull of other bodies, the atmospheric 

^ For an introduction into the Hori-Deprit method see Boccaletti & Pucacco (2002) and Kholshcvnikov 
(1975, 1985). Kinoshita (1977) referred only to the work by Hori (1966). 
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Fig. 1. The perturbed orbit is a set of points belonging to a sequence 
of confocal instantaneous ellipses that are not supposed to be tangent or 
even coplanar to the orbit. As a result, the physical velocity r (tangent to 
the orbit) differs from the Keplerian velocity g (tangent to the ellipse). To 
parameterise the depicted sequence of nonosculating ellipses, and to single it 
out of the other sequences, it is suitable to employ the difference between r and 
g, expressed as a function of time and six (nonosculating) orbital elements: 
^{t , Ci , . . . , Ce) = r(t , Ci , . . . , Ce) — g(t , Ci , . . . , Cq) . In the literature, 
^{t , Ci , . . . , Ce) is called gauge function or gauge velocity or, simply, gauge. 




Fig. 2. The orbit is represented by a sequence of confocal instantaneous 
ellipses that are tangent to the orbit, i.e., osculating. Now, the physi- 
cal velocity r (tangent to the orbit) coincides with the Keplerian veloc- 
ity g (tangent to the ellipse), so that their difference vanishes everywhere: 
$(t , Ci , . . . , Cg) = . This is the so-called Lagrange constraint or Lagrange 
gauge. Orbital elements obeying it are called osculating. 
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and radiation-caused drag, the relativistic corrections, and the non-inertiahty of the reference 
system. 

On Fig. 1 the orbit consists of points, each of which is donated by a representative of 
a certain family of "simple" curves (confocal ellipses). These instantaneous ellipses are not 
supposed to be tangent or even coplanar to the orbit. As a result, the physical velocity r 
(tangent to the orbit) differs from the Keplerian velocity g (tangent to the ellipse). To 
parameterise the depicted sequence of nonosculating ellipses, and to single it out of the other 
sequences, it is suitable to employ the difference between r and g, expressed as a function of the 
time and the orbital elements: $(t , Ci , . . . , Cq) = r{t , Ci , . . . , Cq) — g(t , Ci , . . . , Cq) . 
Evidently, 

i.e., the unperturbed Keplerian velocity is g = dr/dt , while the said difference $ is the 
convective term that emerges when the instantaneous ellipses are being gradually altered by the 
perturbation (and when the orbital elements become time-dependent): $ = ^ [dr/dCj) Cj . 
When one fixes a particular functional dependence of $ upon time and the elements, this 

— * 

function, $(t , Ci , . . . , Cq) , is called gauge function or gauge velocity or, simply, gauge. 

On Fig. 2, the perturbed orbit is represented with a sequence of confocal instantaneous 
ellipses that are tangent to the orbit, i.e., osculating. Under this choice, the physical velocity 
r (tangent to the orbit) will coincide with the Keplerian velocity g (tangent to the ellipse), 

— * 

so that their difference $(t Ci , . . . , Cq) will vanish everywhere: 

$(t, Ci,...,C6)= r{t , Ci, . . . , Cq) - i{t , Cu ■ ■ ■ , Cq) = = 0. 

i=i 

This, so-called Lagrange constraint or Lagrange gauge, is the necessary and sufficient condi- 
tion of osculation of the orbital elements Cj (Brouwer & Clemence 1961). Historically, the 
first attempt of using nonosculating elements dates back to Poincare (1897), though he never 
explored them from the viewpoint of a non-Lagrange constraint choice. (See also Abdullah & 
Albouy (2001), p. 430.) Parameterisation of nonosculation through a non-Lagrange constraint 
was offered in Efroimsky (2002a,b). 

Similarly to orbital dynamics, in attitude dynamics, a complex spin can be presented as a 
sequence of instantaneous configurations borrowed from a family of some "simple rotations" . 
(Efroimsky 2004) It is convenient to employ in this role the motions exhibited by an unde- 
formable free top experiencing no torquesjl Each such undisturbed "simple motion" will be a 
trajectory on the three-dimensional manifold of the Euler angles (Synge & Griffith 1959). For 
the lack of a better term, we shall call these unperturbed motions "Eulerian cones," implying 
that the loci of the rotational axis, which correspond to each such non-perturbed spin state, 
make closed cones (circular, for an axially symmetrical rotator; and elliptic for a triaxial one). 
Then, to implement a perturbed motion, we shall have to go from one Eulerian cone to an- 
other, just as in Fig. 1 and 2 we go from one Keplerian ellipse to another. Hence, similar to 

Here one opportunity will be to utilise in the role of "simple" motions the non-circular Eulerian cones 
described by the actual triaxial top, when it is unforced. Another opportunity will be to use, as "simple" 
motions, the circular Eulerian cones described by a dynamically symmetrical top (and to treat its actual 
triaxiality as another perturbation). The main result of our paper will be invariant under this choice. 
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those pictures, a smooth "walk" over the instantaneous Eulerian cones may be osculating or 
nonosculating. 

The torques, as well as the actual triaxiality of the top and the non-inertial nature of the 
reference frame, will then act as perturbations causing this "walk." Perturbations of the latter 
two types depend not only upon the rotator's orientation but also upon its angular velocity0 



2.2 Delaunay and Andoyer 

In orbital dynamics, we can express the Lagrangian of the reduced two-body problem via the 
spherical coordinates qj = {r , ip , 6 } , then derive their conjugated momenta pj and the 
Hamiltonian 7i(g, p) , and then carry out the Hamilton- Jacobi procedure (Plummer 1918), to 
arrive at the Delaunay variables 

{Qi, Q2, Qs; Pi, P2, P3} = {L, G, H; lo, g, h} = 

(1) 

{ y/JIa , A//ia (1 — e^) , ^/ fia{l — e^) cos i ; — Mo , — , — } , 

where denotes the reduced mass. 

Similarly, in rotational dynamics one can define a state of a spinning top by the three Euler 
angles qj = {ip , 9 , } and their canonical momenta pj = {p^ , Pq , p^} and then carry 
out a canonical transformation to the Andoyer element^ {l,g,h]L,G,H}. By definition, 
the element G is the magnitude of the angular-momentum vector, L is the projection of the 
angular-momentum vector on the principal axis I03 of the body, while H is the projection of 
the angular-momentum vector on the §3 axis of the inertial coordinate system. The variable h 
conjugate to H is the angle from the inertial reference longitude to the ascending node of the 
invariable plane (the one perpendicular to the angular momentum). The variable g conjugate 
to G is the angle from the ascending node of the invariable plane on the reference plane to 
the ascending node of the equator on the invariable plane. Finally, the variable conjugate to 
L is the angle / from the ascending node of the equator on the invariable plane to the the 
hi body axis. Two auxiliary quantities defined through 

T ^ T ^ 

have obvious meaning: / is the angle between the angular-momentum vector and the §3 
space axis, while J is the angle between the angular-momentum vector and the b3 principal 
axis of the body, as depicted on Fig. 3. 

Andoyer (1923) introduced his variables in a manner different from canonical constants: 
while his variables G , H , h are constants (for a free triaxial rotator), the other three, L , I , g , 
do evolve in time, because the Andoyer Hamiltonian of a free top 

When we study the Earth rotation relative to the precessing plane of the Earth orbit about the Sun, the 
frame precession gives birth to a fictitious torque (sometimes called "inertial torque" ) that depends upon the 
Earth's angular velocity. 

^ In attitude dynamics, the Andoyer elements I , g , h play the role of coordinates, while L , G , H are 
their conjugate momenta. In the orbital case, the Delaunay variables L , G , H play the role of coordinates, 
while I , g , h defined as in ([1]) act as momenta. Needless to say, this is merely a matter of convention. (See 
formulae (9.31 - 9.32) in Goldstein 1981.) For example in some textbooks on orbital mechanics the Hamiltonian 
perturbation is deliberately introduced with an opposite sign, while the Delaunay elements I , g , h , too, are 
defined with signs opposite to given in ^ . Under such a convention, the Delaunay elements I , g , h become 
coordinates, while L , G , H act as momenta. 
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S3 




Invariable plane 
(orthogonal to G) 



Fig. 3. A reference coordinate system (inertial or precessing) is constituted 
by axes si , S2 , S3 . A body-fixed frame is defined by the principal axes 
bi , b2 , bs . The third frame is constituted by the angular-momentum vector 
G and a plane orthogonal thereto (the so-called invariable plane). The lines of 
nodes are denoted with i, 1, j . The attitude of the body relative to the reference 
frame is given by the Euler angles hf , If , (j)f . The orientation of the invariable 
plane with respect to the reference frame is determined by the angles h and / . 
The inclination / is equal to the angle that the angular-momentum vector G 
makes with the reference axis S3 . The angle J between the invariable plane and 
the body equator coincides with the angle that G makes with the major-inertia 
axis b3 of the body. The projections of the angular momentum toward the 
reference axis S3 and the body axis b3 are H = G cos / and L = G cos J . 
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2 C 



is a nonvanishing function of the variables /, L and G. (Notations A, B, C stand for the iner- 
tia matrix' principal values that are assumed, without loss of generality, to obey the inequality 
A < B < C .) So, to make our analogy complete, we may carry out one more canonical 
transformation, from the regular Andoyer set {I , g , h , L , G , H} to the "modified Andoyer 
set" { lo , go , h ; Lo , G , H } , where Lq , /q , go are the initial values of L , I , g . The modified 
set consists only of constants of integration, wherefore the appropriate Hamiltonian becomes 
nil§ Therefore, the modified Andoyer set of variables is analogous to the Delaunay set with 
lo = — Mo , while the regular Andoyer elements are analogous to the Delaunay elements with 
I = — M used instead of 4 = — Mo . We would stress that, in analogy with the orbital case, 
the variables h , G , H are constants (and ho = h , Go = G , Ho = H) only in the unperturbed, 
free-spin, case. 

To summarise this section, in both cases we start out with 

. _ m^°^ . _ _ dH^"^ 

^ dp ' ^ dq ^ 

q and p being the coordinates and their conjugated momenta, in the orbital case, or the Euler 
angles and their momenta, in the rotation case. Then we switch, via a canonical transformation 



q = fiQ, P, t) 
p = xiQ, P, t) 



(3) 



to 



0=^=0 . P=-_ = , W=0. (4) 

where Q and P denote the set of Delaunay elements, in the orbital case, or the (modified, 
as explained above) Andoyer set { lo , go , h ; Lo , G , H }, in the case of rigid-body rotation. 

This scheme relies on the fact that, for an unperturbed Keplerian orbit (and, similarly, for 
an undisturbed Eulerian cone), its six-constant parameterisation may be chosen so that: 

1. the parameters are constants and, at the same time, are canonical variables {Q , P} with 
a zero Hamiltonian: 7i*{Q, P) = 0; 

2. for constant Q and P , the transformation equations ([3]) are mathematically equivalent 
to the dynamical equations ([2]). 

In practice, this scheme is implemented via the Hamilton- Jacobi procedure. 



^ We shall not write down the explicit form of this transformation, because it is sufhcient for us to know that 
it is canonical. This follows from the group property of canonical transformations and from the fact that the 
transformations from {l,g,h,L,G,H} to { (p , 9 , ip , Pv ^ PS t Pi' } and from {lo , go , h; Lo , G , H} to 
{l,g,h,L,G,H} are canonical. The latter transformation is canonical, for it is simply the time evolution. 
This canonical transition from Andoyer-type variables to their initial values is not new - see Fukushima & 
Ishizaki (1994). Historically, the first set of rotational elements was constituted by constants (Richelot 1850). 
Serret (1866) found the generating function of a canonical transformation from { ip , 6 , ip , p^^ , pg , } to that 
set. His development was polished by Radau (1869) and Tisserand (1889). The Serret-Richelot set consisted 
of the following constants: { g,, , ft. , — to] G , H , T^in}, where h, G and H coincide with the appropriate 
Andoyer elements, Tj-in is the rotational kinetic energy, to is the initial moment of time, and go is the initial 
value of the Andoyer element g. 
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2.3 Canonical perturbation theory: canonicity versus osculation 

Under perturbation, the "constants" Q, P begin to evolve so that, after their insertion into 

q = /(g(t), P(t), t) , 

(5) 

p = x{Q{t), P{t), t) 

(/ and X being the same functions as in ([3]) ), the resulting motion obeys the disturbed equations 

. _ d + AH) ^ _ d {n^°^ + AH) 

dp ' dq 

We also want our "constants" Q and P to remain canonical and to obey 

• _ d {H* + AH*) ■ _ _ din* + AH*) 

^ ~ dP ' dQ ' ^ ^ 

where 

n* = Q and AH* (Q , P t) = An{q{Q,P,t) , p{Q,P,t) , t) . (8) 

Above all, we wish that the perturbed "constants" Cj = Qi , Q2 , Qs , Pi , P2 , P3 (the De- 
launay elements, in the orbital case, or the modified Andoyer elements, in the rotation case) 
remain osculating. This means that the perturbed velocity will be expressed by the same 
function of Cj{t) and t as the unperturbed one used to. Let us check to what extent this 
optimism is justified. The perturbed velocity reads 

g = g + $ , (9) 

where 

. ,10) 

is the functional expression for the unperturbed velocity; and 

j=l -J 

is the convective term. Since we chose the "constants" Cj to make canonical pairs [Q, P) 
obeying ([7] - [S]) with vanishing Ti* , then insertion of ([7]) into ([TT]) will result in 

^ dQn ^"^^ ^ ^ 5P„ dp ■ ^''> 

n=l n=l 

We see that in some situations the canonicity requirement is incompatible with osculationj^ 
To be specific, under a momentum- dependent perturbation we still can use ansatz ([5]) for 
calculation of the coordinates and momenta, but cannot impose the osculation condition $ = 



^ For the first time, this observation was made in Efroimsky & Goldreich (2003). 
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(i.e., we cannot use g = g for calculating the velocities). Instead, we must use ([9]) with 
the substitution (fT2|) . This generic rule applied both to orbital and rotational motions. Its 
application to the orbital case is illustrated by Fig. 2. There, the constants Cj = {Qn, Pn) 
parameterise instantaneous ellipses which, for nonzero $ , are not tangent to the trajectory. 
In orbital mechanics, the variables preserving canonicity at the cost of osculation are called 
"contact elements" (term coined by Victor Brumberg). The osculating and contact variables 
coincide when the disturbance is velocity-independent. Otherwise, they differ already in the 
first order of the time-dependent perturbation (Efroimsky & Goldreich 2003, 2004). Luckily, 
in some situations, their secular parts differ only in the second order (Efroimsky 2005), a 
fortunate circumstance anticipated by Goldreich (1965), who came across these elements in a 
totally different context unrelated to canonicity. 

The case of rotational motion will parallel the theory of orbits. Now, instead of the instan- 
taneous Keplerian conies, one will deal with instantaneous Eulerian cones (i.e., with the loci 
of the rotational axis, corresponding to non-perturbed spin states). Indeed, the situation of an 
axially symmetric unsupported top at each instant of time is fully defined by the three Euler 
angles g„ = (p , 9 , ip and their time derivatives = (j) , 9 , ip- The evolution of these six 
quantities is governed by three dynamical equations of the second order (the three projections 
of dL/dt = T, where L is the angular momentum and r is the torque) and, therefore, this 
evolution will depend upon the time and the six integration constants: 

Qn = fn {Ci , . . . , Cq , t) , 

(13) 

Qn = g„ (Cl , ■ ■ ■ , Cq , t) , 

where the functions g„ and fn are interconnected via g„ = dfn/dt, for n = ip , 9 , (j). 
Under disturbance, the motion will be altered: 

Qn = fniC,{t), at) , t) , 

(14) 

Qn = gn (Clit) , . . . , Ceit) , t) + (Ci(t) , . . . , at) , t) , 

where 

$.(Ci(t), ...,C6(t),t) = EI^C^. ■ (15) 

j=i J 

If we want the "constants" Cj to constitute canonical pairs {Q, P) obeying ([7]-[H]), then 
insertion of ([7]) into (IT3]) will result in 

so that the canonicity requirement ([7]- [8]) violates the gauge freedom in a non-Lagrange fashion. 

To draw this subsection to a close, let us sum up two facts. First, no matter what the 
Hamiltonian perturbation is to be, the Delaunay (in the orbital case) or the modified Andoyer 
(in the attitude case) variables Q , P always remain canonical. They do so simply because 
they are a priori defined to be canonical - see equations (jl]) and ([7] - [8]) above. Second, as 
we have seen from (fTSl - [T6l) . the osculating character of the Q , P variables is lost under 
momentum-dependent perturbations of the Hamiltonianjf] 



It is possible, of course, to choose the other way and preserve osculation at the cost of canonicity. In the 
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2.4 From the modified Andoyer elements to the regular ones 



So far our description of perturbed spin, ([T3]- [16]), has merely been a particular case of the 
general development (P- [T^ . The sole difference was that the role of canonically-conjugated in- 
tegration constants C = {Q, P) in ([131 -[H]) should be played not by the Delaunay variables 
(as in the orbital case) but by some rotational elements - like, for example, the Richelot- 
Serret variables (see the footnote in subsection 2.2 above) or by the modified Andoyer set 
{lo , Qo , h] Lo , G , H ) consisting of the initial values of the regular Andoyer elements. The 
developments conventionally used in the theory of Earth rotation, as well as in spacecraft at- 
titude engineering, are almost always set out in terms of the regular Andoyer elements, not in 
terms of their initial values (the paper by Fukushima & Ishizaki (1994) being a unique excep- 
tion). Fortunately, all our gadgetry, developed above for the modified Andoyer set, stays appli- 
cable for the regular set. To prove this, let us consider the unperturbed parameterisation of the 
Euler angles g„ = {(j) , 9 , ip) via the regular Andoyer elements Aj = ( / , g , h; L , G , H) : 

qn = fAA.iCt), A,{C,t)) , (17) 

each element Ai being a function of time and of the initial values Cj = {lo , Qo , h; Lo , G , H) . 
When a perturbation gets turned on, the parameterisation ([TTIl stays, while the time evolution 
of the elements Ai changes: beside the standard time-dependence inherent in the free-spin 
Andoyer elements, the perturbed elements acquire an extra time-dependence through the evo- 
lution of their initial values]^ Then the time evolution of an Euler angle g„ = {(^ , 9 , tjj) 
will be given by a sum of two items: (1) the angle's unperturbed dependence upon time and 
time-dependent elements; and (2) an appropriate addition $„ that arises from a perturbation- 
caused alteration of the elements' dependence upon the time: 



The unperturbed part is 




while the convective term is given by 

orbital case, one should simply set $ = in equations (52 - 57) of Efroimsky (2006). In the attitude case, 
though, this will be a more cumbersome construction, never implemented in the literature hitherto. 
^ This is fully analogous to the transition from the unperturbed mean longitude, 

M {t) = Mo + n {t - to) , with Mo, n, to = const , 

to the perturbed one, 

M{t) = Mo{t) + n{t') dt' , with to = const , 

in orbital dynamics. 
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Vf^^ 0- + \-(^JA p.= 9An{q,p) 



where the set Cj is spht into canonical coordinates and momenta hke this: Qj = {lo , Qo , h) 
and Pj = {Lo , G , H) . In the case of free spin they obey the Hamilton equations with a 
vanishing Hamiltonian and, therefore, are all constants. In the case of disturbed spin, their 
evolution is governed by ([7] -[8]), substitution whereof in (120|) will once again take us to (fT6|l . 
This means that the non-osculation-caused convective corrections to the velocities stay the 
same, no matter whether we parameterise the Euler angles through the modified Andoyer 
elements (variable constants) or through the regular Andoyer elements. This invariance will 
become obvious if, once again, we consider the analogy with orbital mechanics: on Fig. 1, the 
correction $ is independent of how we choose to parameterise the nonosculating instanta- 
neous ellipse - through the Delaunay set with or through the one containing M . 



This consideration yields the following consequences: 



(a) Under momentum- dependent perturbations, calculation of the angular velocities via 
the elements must be performed not through the second equation of (fT3l) but through the 
second equation of (IT^ . with (fT6|) substituted therein. The convective term given by (fT6l) is 
nonzero when the perturbation is angular- velocity-dependent. In other words, under such type 
of perturbations, the canonicity condition imposed upon the Richelot-Serret or the Andoyer 
elements is incompatible with osculation. An example of such perturbation shows itself in the 
theory of planetary rotation, when we switch to a coordinate system associated with the orbit 
plane. Precession of this plane makes the frame noninertial, and the appropriate Lagrangian 
perturbation depends upon the planet's angular velocity. The corresponding Hamiltonian per- 
turbation (denoted in the Kinoshita-Souchay theory by E) comes out momentum-dependent. 
In this theory the Andoyer elements are introduced in the precessing frame, and since the 
precession-caused perturbation is momentum- dependent, these elements come out nonoscu- 
lating. For this reason, their substitution into the undisturbed expressions (2.6) and (6.26 - 
6.27) in Kinoshita (1977) will not render the angular velocity relative to the precessing frame 
wherein the elements were introduced. To furnish the angular velocity relative to that frame, 
these expressions must be amended with the appropriate convective terms. 

(b) The above circumstance, instead of being a flaw of the Kinoshita-Souchay theory, 
turns out to be its strong point. It can be shown that Kinoshita's undisturbed expressions 
for the angular velocity via the elements keep rendering the inertial angular velocity, even 
when the elements defined in a precessing frame are plugged therein. Briefly speaking, we first 
introduce the Andoyer elements in an unperturbed setting (inertial frame) and write down the 
expressions, via these elements, for the Euler angles and velocities relative to the inertial frame. 
Then we introduce a momentum-dependent perturbation, i.e., switch to a precessing frame, 
and in that frame we introduce the Andoyer elements. Insertion thereof into the unperturbed 
expressions for the Euler angles and angular velocities gives us the Euler angles relative to the 
precessing frame and (due to the nonosculating nature of the elements) the angular velocity 
relative to the inertial frame, not to the precessing one0 A proof of this fact will be presented 
in Appendix 1.3. 

^'^ This mishap is an example of osculation loss. We introduce the elements in a certain frame (the precessing 
frame of the orbit), plug them into the unperturbed expressions for the Euler angles and for the angular 
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This fact should not be regarded as a disadvantage of the Kinoshita-Souchay theory, because 
in some situations it is the inertial, not the relative, angular velocity that is measured (Schreiber 
et al. 2004, Petrov 2007). Under these circumstances, Kinoshita's expressions for the angular 
velocity should be employed (as long as they are correctly identified as the formulae for the 
inertial angular velocity). 



3 The angular velocity relative to the precessing frame 

In the theory of Earth rotation, three angular velocities emerge: 

uj^ ' = the relative angular velocity, 

i.e., the body's angular velocity relative to a precessing orbital frame; 

/X = the precession rate of the orbital frame with respect to an inertial one; 

uj = the mertial angular velocity, 

i.e., the body's angular velocity with respect to the inertial frame. 

Evidently, the latter is the sum of the two former ones: lj = uj + ^t . 

If some day we develop an experimental technique for measuring the Earth's angular velocity 
relative to the precessing plane of its orbit, we shall have to compare the observations with the 

_^{rel) 

theoretical predictions for the directional angles of this, relative, angular velocity u) 

The Kinoshita (1977) theory was created with intention to furnish the precessing-frame- 
related directional angle^lj of the Earth figure (formulae (2.3) and (6.24 - 6.25) in Kinoshita's 
paper). This theory also provides precessing- frame- related directional angles of the Earth's 
angular-velocity vector (formulae (2.6) and (6.26 - 6.27) in Ibid.). We prove in Appendix 1.3 
below that, contrary to the expectations, the latter expressions render the directional angles 
not of the relative but of the inertial angular velocity uj^ ^ : 



= / + J ( 1 - ^ - ^ ) [cos^ - e cos(2/ + ^)] , (21) 



C C 



and 



J C C 



K = h + —J - ^) - esmi2l + g)] . (22) 



velocities, and here comes the result: while we obtain the correct values of the Euler angles relative to the 
said frame, we do not get the right values for the angular velocity relative to that frame. (Instead, our 
formulae return the values of the angular velocity relative to another, inertial, frame.) This happens because 
the disturbance, associated with a transition to the precessing frame, depends not only upon the Earth's 
orientation but also upon its angular velocity. Or, stated alternatively, because the appropriate Hamiltonian 
variation ATi. depends upon the momenta p canonically conjugated to the Euler angles q : 

op 

Here and hereafter the term "directional angles" will stand for the longitude of the node and the inclination 
of the plane perpendicular to the Earth figure. An analogous meaning is understood for the directional angles 
of the angular velocity. 
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where the an gles / and J are as on Fig. 3, while e is introduced as a measure of triaxiality 
of the rotatoro 

e = [(1/i?) - (lM)]/2 

- il/C) - [(IM) + (1/5)] /2 ^ ^ 

A , i? , C being the principal moments of inertia]^ 

This is a very nontrivial and counterintuitive fact. On introducing the Andoyer variables 
in the precessing frame of orbit, we plug them into the standard expressions for the orientation 
angles and the angular velocity. Doing so, we naturally expect to obtain the orientation and the 
spin rate relative to that precessing frame. We indeed get the body orientation relative to that 
frame, but the rendered angular velocity turns out to be not the one relative to the precessing 
frame wherein the Andoyer elements were introduced. Instead, the standard formulae give us 
the angular velocity relative to some other frame, the inertial one (as if we had used the Andoyer 
variables defined in the inertial frame). This is an interesting (and still underappreciated by 
mathematicians) internal symmetry instilled into the Andoyer construction: we can go through 
a continuum of Andoyer sets (each set introduced in a different precessing frame), but their 
substitution into the standard formulae for the angular velocity will always return the angular 
velocity relative to the inertial frame. 

A proof of this fact begins with a study of the physical meaning of the Andoyer elements 
introduced in a precessing frame (presented in Appendix A. 1.1). Completion of the proof 
demands a sequence of calculations so laborious that we chose to put them into the Appendix 
(see Appendices A. 1.2 - A. 1.3). This entire situation remarkably parallels a similar episode 
from the theory of Delaunay elements in orbital dynamics (see the end of Appendix A. 1.3). 

Now, what if we want to know the angular velocity relative to the precessing frame, i.e., 

{rel) 

uj ? The precessing-frame-related directional angles of this angular velocity will look as 



and 



= /f""' + Ir"' , (24) 



(rei) (inert) (#) /or\ 

hj. = + , (25) 



the extra terms Ir , hr emerging because, as explained above, one has to add the con- 
vective term $ to the unperturbed velocity g , in order to obtain the full velocity q under 
disturbance. Here q stands for the three Eulerian angled g„ = {hj , If , cpf} defining the 
the orientation of the principal axes of the Earth, relative to the precessing frame, so g„ will 

For the Earth, J ~ 10~^ rad, which justifies the common approximation to write all formulae up to the 
first order in J (Kinoshita 1977). The value of the triaxiality parameter is: e = 3.3646441 x 10~^ (Escapa, 
Getino & Ferrandiz 2002). 

The Earth is assumed to be rigid, and its body axes are chosen to diagonalise its inertia matrix. 

Be mindful that in the physics and engineering literature the Euler angles are traditionally denoted 
with {(j) , 9 , ijj) . In the literature on the Earth rotation, very often the inverse convention, (-tp , 9 , cj)) , is 
employed. In the Kinoshita-Souchay theory, these angles are denoted with {h , I , (j)) . The angles defining 
orientation of the Earth's figure are accompanied with the subscript / and are termed as: (/i/ , // , ^/ ) . The 
directional angles of the Earth's angular-velocity vector are equipped with subscript r . For the relative and the 
inertial angular velocities the angles are denoted with {h^ , , (jj^ ) and {h^ , ,0^ ) , 

accordingly. 
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signify time derivatives of Euler angles relative to this precessing frame. The convective terms, 
entering the expressions for g„ = {hf , If , (pf} , can be calculated using formula (fT6|) - see 
the Appendices 3 and 4 below. The ensuing corrections to the Euler angles determining the 
orientation of the instantaneous spin axis will look as 



(s) C • C 

= — TTi — cos/ sin(/i — Hi) + Hi — [sinvri cos/ cos(/i — Hi) + costti sin/ — sin/] 
L L 

+ 0{J^) + OiM/cu) + 0(($/c^)2) 

and 

C = - vr, ^ cos(^.-nO sinvr, sinjh-U,) ^ ^ 0(($/.r) + 0(J$/.' 

L sm/ L sm/ 

These corrections depend upon two angles that define the orientation of the precessing orbit 
with respect to an inertial frame - the inclination tti and the node Hi . 

Let us make rough numerical estimates for the case of the rigid Earth. Putting together 
(|2T|) . fl2^ . and (|26|1 . we see that in the resulting expression for fj"'^ 

(c c \ c 
1 - — - — j [cosg - e cos (2 / + g)] - ^Vi — cos/ sm{h - Hi) 

C 

+ ni — [sinvTi cos/ cos(/i-ni) + costti sin/ - sin/] + O(J^) + 0(J$/u;) + 0(($/cj)^) 

we have the leading term, / , and three additions - of order J ~ 10^^ , of order Je ~ 10~^ , 
and of ordei0 ^/uj ~ tci/co ~ 10~^. We see that the nonosculation- caused convective terms 
$ provide an effect on the spin-axis orientation, which is of the same order as the Je term 
stemming from triaxiality^^ As Irad ~ 0.2 x 10^ " and J ~ 10~^ , then the J term brings 
into Ir a contribution of an arcsecond order, while the $ and Je terms give corrections of 
order milliacrseconds. We also see that the terms of order J ^/u and those of ($/ct;)^ are 
much less than one percent of a microarcsecond and may be neglected. 
Numerical estimates for the expression 

(re!) 

K = h + —J [l-^-^) [sm, - e sm(2/ + ,) ] 

_ c cosjh - n,) _ n, £ ^^""^ sin(/.-n,) ^ _^ og^/^y) + o( j$/^) 

L sin/ L sin/ 

will be similar. 

Formulae fl26l - [27|l constitute the main result of this paper0 In Appendices 2 - 4 we 
present their derivation based on formulae f|T8|) and fl20|) . It would be important to note 



This estimate ensues from the trivial observation that C / L k, uj ^ . Regarding the numbers: according 
to Lieske et ah (1977) and Seidelmann (1992), tti AT /century, while Hi -870'" /century « -2.4 x 
10~^deg/yr. On the other hand, ui ~ 360 deg / day ^ 1.3 x lO^deg/yr whence ^/u ^ tti/uj ^ 10~^. (We could 
as well have used the lERS value of w « 7.3 x 10~^ rad/s ^ 1.3 x lO"^ deg / century ~ 4.7 x 10^" ''/century.) 

The nutational spectra of these two contributions are, however, quite different (secular vs. periodic). 

Through the medium of equations (|62l- [M)) it is also possible to express these corrections via the Euler 
set, instead of the Andoyer variables. 
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that the resulting corrections acquire the form ([25] - [27]) provided the coordinate system co- 
precessing with the orbit is chosen as in Kinoshita (1977), i.e., by three consecutive Euler 
rotations — Hi) R7v(vri) ^zij^i) , letter Z standing for an inertial axis orthogonal to the 
ecliptic of epoch, N denoting the line of nodes, and 3 being a precessing axis perpendicular 
to the ecliptic of date - see Appendix A. 2. 2. 3. Under an alternative choice of axes within the 
co-precessing frame, expressions for ll''"''' and /i^*^"'' will look differently. For example, a 
transition carried out by only two Euler rotations, RAr(7ri) R^(ni) , as in Appendix A. 2. 2. 2, 
will yield expression f]TTD instead of ([25]), and fITM]) instead of fl27]) . 

In principle, ([25] - 1271) might as well be derived by purely geometrical means, i.e., from the 
formula uj = uj + /X . We however chose the method based on (ITSl) and ( [20]) . because 
this method is fundamental and applicable to any kind of momentum-dependent perturbations 
of the Hamiltonian - for example, to the perturbations caused by deviations from rigidity, as 
studied by Getino & Ferrandiz (1994). A similar situation will emerge in the (yet to be built) 
relativistic theory of the Earth rotation. 



4 Conclusions 

In this article we explained that the unperturbed spin states ("Eulerian cones") play in the 
attitude dynamics the same role as the unperturbed two-body orbits ( "Keplerian conies" ) play 
in the orbital mechanics. Just as the orbital elements parameterising Keplerian conies, the 
rotational elements parameterising Eulerian cones may be either osculating or nonosculating. If 
the perturbation depends upon the velocity (in the orbital case) or upon the angular velocity (in 
the attitude case), the condition of osculation is incompatible with the condition of canonicity. 
In these situations the standard equations furnish the Delaunay (in the orbital case) or Andoyer 
(in the attitude case) elements, which are not osculating, - circumstance important when the 
elements are employed for calculation of the velocity or angular velocity. The functional form 
of the expression for a velocity or an angular velocity through elements depends upon whether 
these elements are osculating or not. 

A remarkable peculiarity is shared by the Delaunay and Andoyer elements. Suppose the 
perturbation is caused by a transition to a precessing frame of reference, and the elements are 
introduced in this noninertial frame. Their substitution into the unperturbed expressions for 
the Cartesian coordinates (or the Euler angles) will render the right position (or the attitude) 
relative to the precessing frame wherein these elements were defined. Now, suppose that we 
impose on our elements the condition of canonicity. Since the frame-precession- caused pertur- 
bation is momentum-dependent, the canonicity condition is incompatible with the osculation 
one. Hence, when our elements are inserted into the unperturbed expressions for the velocity 
or angular velocity, they will NOT return the velocity with respect to the precessing frame. It 
turns out, though, that they will render the velocity relative to the inertial frame. While for 
the orbital case this was proven in Efroimsky & Goldreich (2003, 2004), in the current paper 
we proved this fact also for the attitude case. 

This has ramifications for the Kinoshita-Souchay theory of the Earth rotation. In this 
theory, the Andoyer elements are defined in a precessing frame of the Earth orbit. In Kinoshita 
(1977) these elements were ab initio canonical - simply because Kinoshita obtained them via 
a canonical transformation (see section 3 of his work). As demonstrated in sections 2.3 - 2.4 
of our paper, the by-default-imposed canonicity condition made the elements nonosculating. 
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Insertion of such elements into the unperturbed equations for the angular velocity (formulae 
(2.6) and (6.26 - 6.27) in Kinoshita 1977) does not yield the angular velocity relative to the 
frame wherein the elements were defined (the precessing frame). Rather, the equations will still 
furnish the angular velocity relative to the inertial frame of reference. This way of osculation 
loss might be a flaw of the Kinoshita-Souchay theory, had we expected it to render the angular 
velocity with respect to a precessing frame. In reality, the osculation loss is an advantage of the 
theory, because the presently available experimental technique (Schreiber et al. 2004, Petrov 
2007) provides for the measurement of the angular velocity relative to the inertial frame - the 
velocity furnished by the Kinoshita-Souchay theory. 

In the final section we provide expressions ( 126) - 1271) for the body-frame-related directional 
angles of the planet's angular velocity relative to a frame coprecessing with the planet's orbit. 
The method wherewith we calculate these angles is general and applicable to to any kind of 
momentum-dependent perturbations of the Hamiltonian - for example, to the perturbations 
caused by deviations from rigidity. 
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Appendix 1. 

The Andoyer variables introduced in a precessing frame 



A 1.1 Formalism 

Let us consider an unsupported rigid body whose spin is to be studied in a coordinate 
system, which itself is precessing relative to some inertial frame. The said system is assumed 
to precess at a rate p, so the kinetic energy of rotation, in the inertial frame, is given by 

Tk^n = " 2 V + Aij I (^^ + = - 2^ Ji (^u;^ + /ii j 

i=x,y,z 

(30) 

= 2 ^ (^^"'^ + /^x) + 2 ^ (^r'' + ^^y) (^^"'' + /^^) > 

where = (A , B , C) are the principal values of the inertia matrix of the body, uj^' ' is the 

[rel) 

inertial angular velocity (i.e., the one with respect to an inertial frame), uj is the relative 
angular velocity (i.e., the one with respect to a precessing coordinate system), while p, is the 
rotation rate of the precessing frame with respect to the inertial frame. In ( l30l) . both uj 's and 
fl are resolved into their components along the principal axes bi = x , b2 = y , = z of 
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the rotating body. Expression (!30|) is fundamental and stays, no matter whether ft depends 
on the rotator's orientation, or whether it carries a direct time dependence. 
The role of canonical coordinates will be played the Euler angled 

qn = {hf,If,(pf) (31) 

that map the precessing coordinate basis into the principal body basis. To compute their 
conjugate momenta, let us assume that noninertiality of the precessing coordinate system is the 
only angular-velocity-dependent perturbation. Then the momenta are simply the derivatives 
of the kinetic energy. With aid of the formulae for the body-frame components of the relative 
angular velocity|lj 

(reZ) 

LU^ = hf sin If sin0j + If cos0/ , (32) 
ujy"''^ = hf sin If cos (j)f — if sin (pf , (33) 
u^"'^'' = lif cos If + 4>f , (34) 

we obtain: 

= = ^ {^T^ + /^^) sin// sin0/ + B (^^^"" + /i^) sin// cos0/ + C (^a;^'' + /^^) cos// 

Pij = = A l^uj^ + ^ix) COS0/ - B i^ujy + ^lyj sm0/ . 

P,, = ^ = C (c. + , (3 

These formulae enable one to express the angular-velocity components Ui and the derivatives 
Qn = {hf ^ y 'Pf) via the momenta Pn = {Pf^^ , Pj^ , p^^ ) ■ Insertion of ([3S]-[n7]) into 

3 

^ = 5Z ~ ^ = ^/^h/ + ^fP'f + ^fP^j - T + V{hf, If,^f] t) (38) 

n=l 



We would once again remind that the Euler angles, though normally termed ( , 6* , "0 ) i in the astronomical 
literature are often denoted as {ip , 9 , (p) . In the Kinoshita-Souchay theory notations {h , I , </)) are employed. 
The angles defining orientation of the Earth's figure and of the Earth's angular- velocity vector are accompanied 
with the subscripts / and r, correspondingly: ( /i/ , // , (/)/ ) and (hr , Ir , (jir) ■ The directional angles of the 
inertial angular velocity will be denoted with ( /i^.'''^'^*' , ^ ^(ii.=i-t) ^ ^ ^j^Qgg Qf ^j^g relative velocity will be 

called {h!'^"^^ , I^'"'' , ) . 

It should be emphasised, that the components of the angular velocity a; are related to the body 
axes, but the angular velocity itself is the relative one (i.e., that with respect to the precessing coordinate 
system). Our formulae (|32l- [34)) are analogous to equations (2.4) in Kinoshita (1977). At the initial step of his 

(inert) 

development, Kinoshita used his equations (2.4) to express the inertial angular velocity (what we call uj ) 
via the Euler angles introduced in an inertial frame. Then, on having introduced a precessing orbital frame, 
Kinoshita employed these equations for expressing the angular velocity through the Euler angles introduced 
in a precessing frame. Kinoshita did not explore whether this operation would furnish the relative angular 
velocity (what we call uj ) or still the inertial one. 
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results, after some tedious algebra, in 

n = T + AH 



(39) 



the perturbation A7i consisting of a potential term V (presumed to depend only upon the 
time and the angular coordinates, not upon the momenta) and a precession-generated inertial 
term E : 



AH = V{hf , If , ; t) + E 



where 



— Va. cos if + T). COS( 

sin// V^'^/ V ^'f 



(40) 



COS 0/ 



sm 



J- [vn^ - cos//) - p,^ sin, 



expression equivalent to formulae (24 - 25) in Giacaglia & Jefferys (1971)1^ Now let us employ 
the machinery set out in subsection 2.2. The fact that the Andoyer elements are introduced 
in a noninertial frame is accounted for by the emergence of the //-terms in the expression ( !38l) 
for the disturbance A7i . Insertion of (HQ]) into f l20l) entails: 



(In = gn 



dAH 

dpn 



(41) 



where g„ 



hf , If , (f)f , and the convective terms are given by 
dATt fix sin (j)f + fiy cos0/ 



sin/ 



/ 



(42) 



dAH 



— fix COS (f)f + fly sincpf 



(43) 



dAH 

dp, 



[fix SlYKpf + fly COS (f>f) cot If — flz , 



(44) 



It can be shown (Gurfil, Elipe, Tangren & Efroimsky 2007) that the body-frame-related components gi 
of the angular momentum are connected with the Euler angles and their conjugate momenta through 

sini/) . , , 

Qi — p. + p. cosffl — p, smffl cot 7 

^ sin/ ^'f ^ ^'*f ^ 

whence it can be seen that 



92 = Pu^. - Pi^ sm (j) - p^^ cos (j) cot / 

,y^^^^^^^ .V v^^ii u^^^u .|40| is merely another form of the relation AH = — p. ■ G . The latter can also be 

expressed via the Andoyer variables: 
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AH 



fl G 



/ii \J — L?- sin / — [i2 — L?- cos / 
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where g„ stand for hf , If , (j)f , and Pn signify the corresponding momenta, while fix , fiy , fiz 
are the components of /x in the principal axes of the body. 

A 1.2 The physical interpretation of the Andoyer variables 
defined in a precessing frame 

The physical content of the Andoyer construction built in an inertial frame is transparent: 
see Fig. 3 and explanation thereto. Will all the Andoyer variables and the auxiliary angles 
/ and J retain the same physical meaning if we re-introduce the Andoyer construction in 
a noninertial frame? The answer is affirmative, because a transition to a noninertial frame 
is no different from any other perturbation: precession of the fiducial frame (§1 , §2 , §3) is 
equivalent to emergence of an extra perturbing torque, one generated by the inertial forces 
(i.e., by the fictitious forces emerging in the noninertial frame of references). In the original 
Andoyer construction assembled in an inertial space, the invariable plane was orthogonal to 
the instantaneous direction of the angular-momentum vector: if the perturbing torques were 
to instantaneously vanish, the angular-momentum vector (and the invariable plane orthogonal 
thereto) would freeze in their positions relative to the fiducial axes (§1 , §2 , §3) (which were 
inertial and therefore indifferent to vanishing of the perturbation). Now, that the Andoyer 
construction is built in a precessing frame, the fiducial plane is no longer inertial. Never- 
theless if the inertial torques were to instantaneously vanish, then the invariable plane would 
still freeze relative to the fiducial plane (because the fiducial plane would seize its precession). 
Therefore, all the variables retain their initial meaning. In particular, the variables / and 
J defined as above will be the angles that the angular-momentum makes, correspondingly, 
with the precessing §3 space axis and with the 63 principal axis of the bodyo Among 
other things, this explains why Laskar & Robutel (1993) and Touma & Wisdom (1993, 1994), 
who explored the history of the Martian obliquity, arrived to very close results. Both groups 
rightly used the angle / as an approximation for the obliquity. While Touma & Wisdom 
(1993, 1994) employed (a somewhat simplified version of) the Kinoshita formalism in an in- 
ertial frame, Laskar & Robutel (1993) used this machinery in a precessing frame of the orbit. 
Now we understand why they obtained so close results, with minor differences stemming, most 
likely, from averaging- caused error accumulation in the latter paper. (The computation by 
Touma and Wisdom was based on unaveraged equations of motion, while Laskar and Robutel 
employed orbit-averaged equations.) 

A 1.3 Calculation of the angular velocities via the Andoyer 
variables introduced in a precessing frame of reference 

Let us now have a look at the well-known expressions 

(reZ) 

= hf sin// sin0/ -|- // cos0/ , (45) 

{rel) 

ujy = hf COS0/ sin If — If sin0j , (46) 

(rel) 

= (j)f -\- hf cos If , (47) 

On the interrelation between the Andoyer variables, referred to an inertial frame, and those referred to a 
moving frame see equation (3.3) in Kinoshita (1977). 
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,(rei) 



for the principal-axes components of the precessing- frame-related angular velocity uj . These 
formulae render this angular velocity as a function of the rates of Euler angle's evolution, so one 
can symbolically denote the functional dependence (jl5]-|17]) as uj = uj{q) . This dependence 
is linear, so insertion of (14T!) therein will yield: 



(4J 



with A denoting the set of Andoyer variables, and p signifying the canonical momenta 
corresponding to the Euler angles. Direct substitution of (H^l - IHl) into (H5] - H7|) will then show 
that the second term on the right-hand side in fHHj) is exactly — /x : 



C3{q{A)) = L5(g(A)) - il 



(49) 



Since the is, ab initio, the relative angular velocity uj ' (i.e., that of the body frame 

relative to the precessing frame), and since is the precession rate of that frame with respect 
to the inertial one, then uj{ g{A) ) will always return the inertial angular velocity of the body, 
l3 (i.e., the angular velocity relative to the inertial frame). It will do so even despite the 
fact that now the Andoyer parameterisation is introduced in a precessing coordinate frame! 
In brief, the above line of reasoning may be summarised as: 

QiqiA)) = '^(g(^)) + i^^idAH/dp) , 



Q{qiA)) = 



Gj{dAn/dp) = p: 



__^{rel) __^{inert) 



where the entities are defined as follows: 



^1 I AW _(inert) 



(50) 



the relative angular velocity, 

i.e., the body's angular velocity relative to a precessing orbital frame; 



fl = the precession rate of that frame with respect to an inertial one; 

u) = the mertial angular velocity, 

i.e., the body's angular velocity with respect to the inertial frame. 

This development parallels a situation in orbital dynamics. There the role of canonical 
elements is played by the Delaunay set C = (Q; P) = {L , G , H ; — Mo , —uj, —fl) with 

L = /^i/2«^/2 ^ Q = fi'/^a'/' (1 - e'f' , H = fi'/' a'/^ (l - e^f^' cos « , 

the parameters a, e, i, u, Q, Mo being the Kepler orbital elements. In the unperturbed setting 
(the two-body problem in inertial axes), the Cartesian coordinates r = {xi,X2, Xa) and 
velocities {xi , X2, X3) are expressed via the time and the Delaunay constants by means of 
the following functional dependencies: 

— * — * 

^ = f{C, t) and V = g(C, t) , where g = df/dt . (51) 



21 



If we want to describe a satellite orbiting a precessing oblate planet, we may fix our reference 
frame on the precessing equator of date. Then the two-body problem will get amended with 
two disturbances. One, AHoUate , caused by the presence of the equatorial bulge of the planet, 
will depend only upon the satellite's position. Another one, ATiprecess , will stem from the 
noninertial nature of our frame and, thus, will give birth to velocity-dependent inertial forces. 
Under these perturbations, the Delaunay constants (now introduced in the precessing frame) 
will become canonical variables evolving in time. As explained in subsection 2.3, the velocity- 
dependence of one of the perturbations involved will make the Delaunay variables nonosculating 

— * 

(provided that we keep them canonical). On the one hand, the expression r = f {C{t) , t) 
will return the correct Cartesian coordinates of the satellite in the precessing equatorial frame, 
i.e., in the frame wherein the Delaunay variables were introduced. On the other hand, the ex- 
pression g ( C , t ) will no longer return the correct velocities in that frame. Indeed, according 
to (P- [TO|) . the Cartesian components of the velocity in the precessing equatorial frame will be 
given by g ( C , t) + dAHprecess/dp . However, since the second term of this sum is equal to 
— p, X r, then g ( C , t) turns out to always render the velocity with respect to the inertial 
frame of reference (Efroimsky & Goldreich 2004, Efroimsky 2005). 



Appendix 2. 

The instantaneous angular velocity 
in the Kinoshita-Souchay theory 



The main burden of subsection 2.3 in the text above was to highlight the need to add 
the convective term $ to the unperturbed velocity g , in order to obtain the full veloc- 
ity q under disturbance. Here q stands for a vector consisting of the three Eulerian 
angles g„ = hj , If , (pf defining the orientation of the principal axes of the Earth relative 
to the precessing frame. The corresponding convective terms, entering the expressions for 
qn = hf , If , (pf , are given by formula (l20l) . Our eventual goal will be to calculate the cor- 
responding corrections to the Euler angles determining the instantaneous axis of rotation in a 
precessing frame of reference. 

A 2.1 The unperturbed velocities 

In this subsection we shall write the unperturbed Euler angles' partial time derivatives 
g^ = dqn/dt as functions of these angles and of the Andoyer variables. 

In the Kinoshita-Souchay theory of Earth rotation, the Euler angles defining the figure of 
the Earth are denoted with 

qn = (hf , If , (pf) , (52) 

the subscript standing for "figure." 

Now, let us denote the principal body axes with 1,2,3 and the appropriate moments 
of inertia with A, B, C (so that A < B < C ). The angular momentum L is 
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connected with the Earth- figure Euler angles via the body- frame components (145) - Wf} of the 
inertial-frame-relatecS angular velocity uJ* ' : 



A u. 



(inert) 



A ihf sin If sin (j)f + If cos( 



(53) 



B uj. 



(inert) 

y 

{inert) 
Z 



B ihf sin If cos (j)f — If sini 



C [hf cos// + 



(54) 



(55) 



On the other hand, the body-frame components of the angular momentum will be related to 
the Andoyer elements through!^ 



VG^ _ ^2 sin/ 



(56) 



cos , 



(57) 



L . 



(5^ 



Substituting (1561-I58!) into ([53] -[55]) and solving for the rates of change of the Euler angles will 
entail: 



dh 



f 



1 



dt sin I 



f L 



— sm(j)f + — cos< 



1 



sin // 



_ ^2 



sin I sin (j)f ^ cos/ cos < 



A 



B 



dt 



A 



COS( 



B 



sm I 



VG2 - L2 



sin / cos(/)/ cos / sine/)/ 



A 



B 



dt 



C 



cot If 



'a 



smi 



+ COSi 

B 



L 
C 



L2 cot / 



/ 



sin I sin (J) f ^ cos/ cos < 



A 



B 



where we deliberately replaced hf , If, (f)f with dhf/dt , dif/dt , d(f)f/dt , because so 
far we have been considering the situation of no disturbances turned on (i.e., the case when 



Be mindful that in formulae ([45l - 147 | the notations hf , If , (f)f stood for the Euler angles defining the 
body's orientation relative to a precessing frame. For this reason, ([45l - 147]) furnished the relative angular 

(rei) 

velocity a3 . In this subsection we are beginning with the unperturbed situation, when the orbit frame is 
yet assumed to be inertial. Hence, at this moment, hf , If , (j>f yet denote the angles relative to the inertial 

_^(iTiert) 

frame, and hence the same formulae render the inertial angular velocity uj 

At this point, we are discussing the unperturbed case, with no frame precession. However, as explained 
in subsection Al.2, the interconnection (ISl 



between the Andoyer elements and the components of L 
will stay valid also when the precession is "turned on" (and both the elements and the components of L are 
introduced in a precessing frame of reference). 
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the full derivatives coincide with the partial ones, and lack convective terms). Our next step 
will be to turn on the disturbance A7i , which will include a transition from an inertial frame 
to the precessing frame of the Earth's orbit. In accordance with formulae f|T^ and (fT6!) . this 
transition will generate additions to the derivatives (jSni - EID, the additions that make the 
difference between a total and a partial derivative. 



A 2.2 Turning on the perturbation — switching to a precessing frame 

Our goal here is to derive the convective terms = ( , , ) that are to be added 
to the partial derivatives fl59]- [6T]). to get the full time derivatives g„ = ( /i/ , // , 0/ ). 

A 2.2.1 Generalities 

As explained by Kinoshita (1977), the undisturbed dependence of the Euler angles of the 
Earth's figure upon the Andoyer elements can be approximated with 

hf = h + ^-j sing + 0{J^) , (62) 

= I + J cosg + 0{J'^) , (63) 

(f)f = I + g - J cot I sing + 0{J^) , (64) 

J and / being the angles that the invariable plane (the one orthogonal to the angular momen- 
tum G ) makes with the body equator and with the ecliptic plane of date, correspondingly. 
(For the Earth, J is of order 10^^ . ) As evident from Fig. 3, these angles are interconnected 
with the Andoyer variables L and G through formulae 



G cos J (65) 



H = G cos/ . (66) 



and 



Under perturbations, formulae fl62] - [64l) will stay valid. However, the expressions for the angles' 
evolution rate, ([59]- [61]), will acquire convective additions f l20]) caused by the loss of osculation. 
These additions, entering the expressions for g„ = {hf , If , </>/), will read, accordingly, as 



dAH ^ dAH ^ dAH 

^ ophf ' dpi ^ dp^ 



So our next step will be to calculate these three terms. 

Among the perturbations entering the Kinoshita theory, there is a so-called "E term." It 
emerges due to a transition from an inertial frame to a noninertial one, i.e., from a coordinate 
system associated with the ecliptic of epoch to the one associated with the ecliptic of date. 
Simply speaking, in the Kinoshita theory the Earth rotation is considered in a noninertial 
frame of the terrestrial orbit precessing about the Sun. In Kinoshita (1977), the xy plane of 



24 



this noninertial frame is referred to as the moving plane. In his theory, this "E term" is the 
only one dependent not only upon the instantaneous orientation but also upon the angular 
velocity of the Earth (or, in the Hamiltonian formulation, upon the momenta conjugate to 
the Euler angles of the Earth's figure). Hence, in this situation dATi/dpj = dE/dpj . The 
expressions for A7i and the "E term" are rendered by formulae (15^ - HU]) where p^^ , p^^ , p^^ 
denote the canonical momenta, while Hx , iJ-y , fJ-z signify the body-frame components of the 
angular rate at which the orbit plane is precessing relative to an inertial coordinate systemic 
In order to continue, we need the expressions for the body-frame components fix , fJ^y , fJ-z ■ 
These can be obtained from the precessing-frame components /xi , /i2 , A^s by means of the 
appropriate rotation matrix: 

l^z 

COS0J COS hf — sin (j)f cos If sinhj cos0j sin hf + sin0j cos If cos hf sin0j sin If 

— sin(^j coshf — cos0/ cos// sinhf —sincpf sinhf + cos0/ cos// coshf cos0/ sin/j 

sin I f sink f —sin If coshf cos// 

Since no other contributions in ATi. other than E depend upon the momenta, then 

^ ^ — , . SinCpf COS0f r • r r 7 

<Phf = — = = Lix — Uv = Ui cot If Sinhf — U2 cot if coshf — 

1 a„ sin// sin// 

$7^ = — = — — = — Hx COS0/ + Hy sin0/ = — jii coshf — /i2 sinhf , 



dAH 


dE 


dp., 


dPnj 


dAH 


dE 


% 


dp., 


dAH 


dE 


9P,, 


9P,, 



. sin0/ cos// COS0/ cos// sinhf coshf 

'^</'/ = ^ — = ^ — = '■ — } '■ — F f^y ~ = — ■ — r ^1 — '■ — F' 

J smi/ smi/ smi/ smi/ 

Naturally, none of the $ terms bears dependence upon 0^ . 

To write down the $ terms as functions of the longitude Hi and inclination tti of the 
ecliptic of date on that of epoch, we shall insert into ([69] - [71]) the appropriate expressions 
for /ii , /i2 , fi3 ■ However, at this point care is neededj_because of the freedom of choice of a 
coordinate system co-precessing with the orbital planel 



A 2.2.2 The precession rate p, as seen in a certain coordinate system 
associated with the precessing equator of date 

Let the inertial axes {X , Y , Z) be fixed in space so that X and Y belong to the 
ecliptic of epoch. A rotation within the ecliptic-of-epoch plane by longitude Hi , from the axis 

We would point out that in Kinoshita's theory the origin of both Hi and h is the mean equinox, whereas 
in our formaUsm the origin simply coincides with the x axis. For our Hi and h to coincide with those of 
Kinoshita, not only must we choose our inertial coordinate system with its xy plane being within the ecliptic 
of epoch, but we should also choose the x axis to coincide with the mean equinox of epoch. Similarly, not 
only should our precessing frame to be associated with the ecliptic of date, but the precessing x axis whence 
we reckon the angles should be placed exactly at the angular distance of — Hi from the node - see Fig. 2 in 
Kinoshita (1977). (Mind that in the presence of precession Kinoshita employs notation h' instead of h .) 

We are grateful to Hiroshi Kinoshita who explained to us the choice accepted in his works. 
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X , will define the line of nodes. A rotation about this line by an inclination angle tti will give 
us the ecliptic of date. The line of nodes, 1 , along with axis 2 naturally chosen within the 
ecliptic-of-date plane, and with axis 3 orthogonal to this plane, will constitute the precessing 
coordinate system, with the appropriate basis denoted by ( §1 , §2 , §3 ) . For example, the 
unit vector §3 reads in the inertial axes ( X , y , Z) as 

T 

§3 = ( sin TTi sin Hi , — sin tti cos Hi , cos tti ) . (72) 
The Earth's angular velocity relative to the inertial and precessing axes obey 



-^(inert) -^(rel) 



(73) 



/7 being the precession rate of the precessing axes relative to the inertial axes ( X , Y , Z) . 
In the inertial axes, this rate is given by 



a' 



TTi COS Hi , TTi sin Hi , Hi 



(74) 



because this expression satisfies the equality /x' x §3 = §3 , as can be easily seen from (!72l) 
and dZH). 

In a frame precessing with the ecliptic, the precession rate will be represented by the vector 



where 



R 



e— ►d 



Ri(7ri)Rz(ni) 



COS Hi 



— cosTTi sin Hi 
sin TTi sin Hi 



sin Hi 

cos TTi COS Hi 
— siUTTi COS Hi 



(75) 





sin TTi 

COSTTi 



(76) 



is the matrix of rotation from the ecliptic of epoch to that of date. From f[75] - [76]) we get the 
components of the precession rate@ as seen in the co-precessing coordinate frame (1 , 2 , 3) : 



fj. = [ fll , fl2 , fJ'3 ) = ( TTi , III sinVTi , ni COSTTi 

Substitution of these components into ([69] - [H]) entails: 

= cot// sin hf — Hi sinvri cot// cos/i/ — Hi costti 



(77) 



(78) 



TTi COS hf — 111 sinvTi sink 



f ' 



(79) 



sin /if cos hf ■ 

TTi + — — Hi smvTi . 



sin / 



/ 



sin/ 



/ 



(80) 



Equivalently, one can find the components of /7 as the elements of the skew-symmetric matrix Re^d R-e^d ■ 
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A 2.2.3 The precession rate ft as seen in a different coordinate system 
associated with the precessing equator of date 
(the system used by Kinoshita 1977) 

In the preceding subsection the transition from the ecliptic of epoch to the one of date 
was implemented by two Euler rotations: Re^d = RAr(7ri) ^zij^i) = Ri(7ri) Rz(ni) . 
The axis 1 of the precessing frame was assumed to coincide with the line of nodes, . 
Evidently, this choice was just one out of an infinite multitude. An alternative option was 
employed by Kinoshita (1977), who used a sequence of three Eulerian rotations: R-el^j^ = 
R3( — Hi) R7v(vri) R^(ni) . Specifically, having performed the two rotations described above, 
Kinoshita then rotated the axis 1 within the ecliptic of date by angle — Hi away from the 
line of nodes N . Due to reasoning analogous to what was presented in the subsection above, 
the sequence of three rotations gives, instead of (177j) . the following expression^ 

T 

fl = ( ^1 , /i2 , Us) 

T 

= ^TTi COS Hi — n sin Hi sinvri , tti sin Hi + 111 cos Hi sinvri , Hi costti — 111 j 

Insertion thereof into (1691- ITTI) will yield: 

^hf = cot// sm{hf — Hi) — Hi sinvri cot// cos(/i/ — Hi) — Hi costti + Hi , (82) 

= —771 cos(/i/ — Hi) — Hi sinvTi sm(hf — Hi) , (83) 



sm{hf - Hi) . cos{hf - Hi) • . 

$0,- = -TTi + ^-^^ -Ui smvTi . 84 

sm I / sm I / 



(81) 



A 2.3 The perturbed velocities 



According to f|T8|) . to get the full evolution of the figure-axis Euler angles under pertur- 
bation, one should sum the unperturbed velocities, given by the partial derivatives fl5^- IM]) . 
with the appropriate convective terms (1^- ITTj) : 



dt 



(85) 



c 



sin / 



VG2 



L2 



/ 



sin / sin (/)/ ^ cos / cos ( 



A 



B 



+ TTi cot// sin/i/ — Hi siuTTi cot// cos/i/ — Hi cos 



TTl 



In the precessing frame, the angular momentum reads: (G sin/ sin/i, — G sin/ cos/i, H) , quantities 
I , h , and H being as in Fig. 3. This, together with ([51]) and the formula AH = E = p, ■ G , will entail: 

AH = E = ill H {1 - cosvTi) - TTi G sin/ sm{h - Hi) + 111 G sin/ cos(/i - Hi) sIutti , 

which coincides with expression (3.4) in Kinoshita (1977). 
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'dt 



c 



VG2 - L2 



sin / cos ( 



cos / sin I 



TTi COS hf — 111 sinvTi sin hf 



16) 



(j)f 



dt 



+ $ 



c 



L 

C 



- - L2 cot If 



sin / sin ( 



+ 



cos / COS ( 



B 



sin h f cos /i f • 

TTi + ^ — Hi smTTi 



sin/ 



/ 



sin I 



f 



These expressions will help us to determine the instantaneous orientation of the spin axis. 



^7) 



A 2.4 The precessing-frame-related angular velocity 
expressed through the Andoyer elements 
introduced in the precessing frame. 

Let angles h^^^''' and J,^"^"'' be the precessing-frame-related node and inclination of vector 
of the angular velocity relative to the precessing frame, and let u 
relative spin rate. Ou 
of hf,If,(j)f,hf, if 



(rel) 



be the 

relative spin rate. Our next step will be to derive expressions for h^^"'''' and I^J"''^ as functions 

and then to substitute the expressions ([85]- [87]) for hf , If , (pf 

f , If , (pf only. To accomplish this step, let 



therein, in order to express h^J^"'^ and /^'^"'^ via hf , I 
us begin with the formulae interconnecting the precessing-frame components of the angular 
velocity relative to the precessing frame with the figure-axis Euler angles and with the spin-axis 
Euler angles. These formulae are fundamental and perturbation- invariant: 



CO 



(rel) 



If COS hf + 



f 



sin If sin h 



f 



(88) 



a;. 



(rel) 



If sin h 



f 



sin If cos h 



f 



(89) 



CO. 



(rel) 



hf + (j)f cos// 



and 



(rel) 



{rel) , {rel) , - (rel) 

uj sm L sm h^ 



(90) 



(91) 



UJ. 



(rel) 



{rel) . {rel) - {rel) 

UJ sm L cos h^ 



(92) 



UJ 



(rel) 



{rel) ^{rel) 

UJ cos K 



(93) 



Both the inertial and relative spin rates, and dj''^'''' , can be most conveniently calculated 

in the body frame. In that frame, the inertial angular velocity can be written as 

T T 



{inert) 



Lx 



B ' C 



G sin J sin I 
A 



G sin J cos / 
B 



G cos J 
C 



(94) 
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whence its absolute value turns out to be 



u^'-' = ^l i^X + (^X + (I'' 



To derive the relative rate, square the obvious equality uj = oj — n , to obtain: 

, ^(rel) ,2 ^(inert) , 2 ^(inert) ^ ^2 

{uj ) = {uj ) - 2uj ■ fi + p, . (96) 



Hence, 

= - a + 0[{^/u;f)] = ^[l - a + 0[{^/u;f + 0{J'))] (97) 



and 

C 



^ = + « + O (($/.; )^)] = [1 + « + 0(($/u;)^) + 0(J^)] , (98) 



a; 



where 



a = " ^ . (99) 

(a; ) 



is of order ^/u . Dot-multiplying (1941) by (1681) . we arrive at: 

G cos J liz ^/ t:^ / X 

^ (a; ) 

C 

= — [ni sin// smhf — /i2 sin// cos/i/ + /^a cos//] + 0(J$/c<j) + O(J^) . (100) 
ij 

In the coordinate system described in A. 2. 2. 2, substitution of fl77j) makes the latter read 

a = — (^TTi sin// sin/i/ — Hi sinvri sin// cos/i/ + 111 costti cos//j + 0(J$/u;) + O(J^) . (101) 

In Kinoshita's coordinate system described in subsection A. 2. 2. 3, one should use for the com- 
ponents of fx not expression (1771) but (IHTl) . insertion whereof into (llOOj) results in: 

a = ^ (^^1 sin// sin(/i/ — Hi) — Hi sinvri sin// cos(/i/ — Hi) -|- Hi costti cos// 

(102) 

- Ill cos//) + 0(J$/cu) + O(J^) . 

This formula will enable us to derive approximate (valid up to O {.P) + {J^/uj) + {{^/ujY) ) 
expressions for /i^.'"'' and /^'^'''' expressed as functions of the Andoyer variables. 
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Appendix 3. Expression for / 



(rel) 



Together, (!90!) and (!93|) will give: 
or, equivalently. 



(7'eZ) {rel) 

uj cos L. 



(rel) 



cosL. — cos if 



hf + 



hf + 



COS/f 



— u; I COS-// 



(103) 



(104) 



Since we are planning to carry out all calculations neglecting terms O ( J^) , and since the 
three inclinations // , Ij"^ , I differ from one another by quantities of order 0{J) , we can 
approximate the left-hand side of fll04p with the first-order terms of its Taylor expansion: 



(rel) . . 
UJ — Sm if 



Arel) 



wherefrom 



j(rel) 

K - h 



(rel) 



sin If 



+ 



UJ 



(rel) 



hi 



- 1 



— UJ 



(rel) 



cos// + O {.f) , (105) 



cot //) + O ( J2) 



UJ 



(rel) gjj^ J 



+ 



/ 



dt 



UJ 



(rel) 



- 1 



cot//) + O (J2) 



To get rid of the time derivatives, employ formulae (j85]-[87]) or (1591- [6T!): 



j(rel) 

i^ —if 



(rel) ■ 2 T 

UJ sm L 



cot If 



VG2 - L2 



UJ 



(rel) 



- - VG^ - L2 cot// 



sin / sin i 



+ 



cos / COS ( 



sin / sin ( 



+ 



cos / COS I 



B 



(106) 



(107) 



(rel) 

UJ sm 



— (^f,^ + cos//) + O (J2) 



For (pf we can use the approximation 0/ = / + g — J cot / sing + 0{.P) . Besides, in the 
terms of order J and of order $/a; we can substitute cu*'^'''' with = L/C + 0{J^) . 

Such substitutions will entail errors of orders 0{J(^/uj) and 0( ("I'/a;)^ ) . However, in the 
leading term we must use fl97l- llOll) . Thus we get: 



(rel) 



COt/f 1 



UJ 



1 L \ - L2 , 2 ^ 



c 



UJ 



sin^ // 



sin / sin i 



+ 



cos / COS ( 



B 



1 1 



a;''"'' sm i/ 



— - + cos//) + 0(J2 
m /f ^ ^ ' 



cot /f 1 - 



1 + a L 



UJ 



{inert) 



G sin J 

G/C 



sin / r • \ cos / , , 
sm (/ + g — J cot i smg) H cos ( / 



A 



B 



g — J cot / sin J 
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1 1 
TJc^ 



- cos//) + 0{.f) + 0(J$/^) + 0(($/cj) 

1 f 



— a cot If — C J 



sin / cos / 

sm(/ + g) + cos(£ + 



5 



C 1 



whence, by using f llOip and the formula If = I + J cosg + 0{J^) , we arrive at: 

i f C C 1 

ij" = I + J < cos g ~ ^ sin / sin (/ + g) — — cos / cos (/ + g) > — a cot // 



108) 



- J ($ + <f<^, cos//) + O (J2) + O (J$/a7) + O ( ) 

Via trigonometric transformations, the second term gets simphfied as: 

C C C cosg - cos(2/ + g) 

cosg — - sm / sm (t + gj — — cos / cos (/ + gj = cosg — — 

A Id a 2 



(109) 



C cosg + cos(2 / + g) 
B 2 



C C 

^ ~ ~ 25 ' ~ e cos(2/ + g); 



Insertion of fl^ - [7T]) into the fourth term will make it look: 

- T—r 0^^/ +"^0/ cos// = - 

L sm If ^ ^ ' ' L sm // 



Hence, we have for /, 



(110) 



a cot / 



/ 



+ 0(J2) + 0{J^/u) + 0{{^/uf) 



111) 



where the parameter e , given by (1251) . is the measure of triaxiality of the rotator. 

In the precessing coordinate system obtained from the inertial one by two Euler rotations, 
as in subsection A. 2. 2. 2, we must now substitute (11011) for a and (177|1 for /is , to get: 

(c c \ c 
1 - ^ - ^ j [ COS g - e cos (2 / + g) ] - — TTi cos // sink 

C 

+ — III ( sinvTi cos// cosh + costti sin// ) + O (J^) + 0{J^/uj) + 0{{^/ujf) . 



(112) 
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In the Kinoshita precessing axes obtained from the inertial ones by three rotations, as in 
subsection A. 2. 2. 3, we should substitute (11021) for a and (IHT!) for /is , to get: 

^ ~ 2A ~ 2b) [ g ^ ^ cos (2 ^ + g) ] - TTi cos // sm{h - Hi) 

C 

+ — III (sinvTi cos// cos{h - Hi) + costti sin J/ - sin//) + O (J^) + 0(J$/u;) + 0{ [^/uj^ ] 
ij 

To arrive to the final expression for J^'^"'' , we shall make use of ( 162) - [6^ . These formulae will 
enable us to substitute, in the above expression. If and hf with / and h, correspondingly. 
All in all, in the coordinate system as in A. 2. 2. 2 we have: 



(113) 



/^"'^ = ^ + y- ~ [c°sg ~ e cos (2 / + g) ] - ^ TTi cos/ sin/i 

C 

+ — fli ( sinvTi cos/ cos/i + costti sin / ) +0 (J^) + 0(J$/u;) + 0(($/u;)^) . 

In the coordinate system as in A. 2. 2. 3, we obtain: 

(C C \ C 

^ ~ 2A ~ w) [ cos g ^ ^ cos (2 ^ + g) ] - TTi cos / sm{h - Hi) 

C 

+ — 111 [sinvTi cos / cos(/i - Hi) + cosvti sin / - sin/] + O (J^) + 0(J$/cj) + 0( i^/uj'f\ 

In (I114p and (11151) . the first two terms coincide with those given by the second of formulae 
(2.6) in Kinoshita (1977). They make j^™'"''' ^ while the third term is /^*' . 



Appendix 4. Expression for 



(114) 



;ii5) 



Expressions ( 1881) and ( I9T1) result in 

u;''^'''' sin /^'^'''' sin h!'^"''^ = If cos hf + 4>f sin // sin hf , (US) 

while dHHD and ([92]) entail 

— tu' ' sin /j', ' cos hi ' = If sin hf — 4>f sin // cos /i/ . (117) 

Multiplying the former with cos hf and the latter with sin hf , and then summing up the two 
results, we arrive at 

If = uj smi^ sm h^ — ht ] . 118 
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Since the difference /i^'^"'' — hf is expected to be of order 0(J) + 0{^/uj) , the above formula 
may be rewritten as 



(rei) 

— 111 



( rei ) • T 

uo smij 
or, according to (198!) . as 

1 + a 1 /(9J 



+ 0(J2) + 0{{^/ujf) + 0(J<l>/cu) 



(119) 



(rei) 



u;''"^'''' sin 



+ <f/,) + O(J^) + 0((<|./c.)2) + 0(J$/c.) (120) 



where a is of order 0($/u;) and therefore may be dropped. Recall that, according to fl95|) . 
the absolute value of the angular-velocity vector can be approximated, up to O(J^) , with 
^(»nert) ^ ^ G / C , whllc a/G^ " L"^ cdji bc expressed as ^G^ — U- = G sin J ^ 

G J L J. Together with approximations 0/ = / + g — J cot / sing + O(J^) and // = 
I + J cosg + O(J^) , it will enable us to rewrite (|60l) as 

sin / cos(/ + g) cos / sin(/ + (?) 



dt 



L J 



A 



B 



+ 0(J2 



(121) 



insertion whereof in (11211) will then entail 



— fit 



JC 
sin If 



sin / cos(/ + g) cos / sin(/ + g) 



A 



B 



C 1 



L sin/ 



/ 



+ 0( J2) + 0( {<^/uof ) + 0( J $/cu ) (122) 



This, along with hj = h + J sin (7/ sin / + 0{.P) and If = I + J cosg + 0{.P) , yields: 



(rei) .sing JC 

/i„ = h + J - — - + - — - 
sm I sm I 



sin / 
~A~ 



cos (/ + g) 



COS / 



sin (/ + g) 



+ T + 0( 

L smi 



J / . C sin(2/ + g) - sing 

sm i A 2 



C sin(2/ + g) + sin) 
B 2 



+ (-/ii coshf - /i2 sin/i/)+ 0(J2) + 0(($/cu)') + 0(J$/. 

Ij sm J 



c 

h + 1 — — - — — I I smg 

sin / V 2A 2B ' ^ ^ 



e sin (2 / + g) 



+ 



C 1 
L sin/ 



-/ii coshf - fi2 sin hf)+ 0{J^) + 0((<l>/cu)2) + 0{J<^/lj' 
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To get the answer in the precessing coordinate axes obtained from the inertial ones by two 
Euler rotations, as in A. 2. 2. 2, we substitute (1771) for /ii and /i2 . It yields: 

(rei) 

K =/^ + — ^ (^l-^-^J [smg- esm(2/ + g)] 

-ir.^'-^-U,^':^^^^^ . (124) 

L smi L smi 

To obtain the answer in the Kinoshita precessing coordinate system obtained from the inertial 
one by three Euler rotations, as in A. 2. 2. 3, we should substitute (ISTj) for /ii and 112 ■ This will 
entail: 

. C cos(h — Hi) • C sinvTi sm(h — Hi) ^, ^o, ^, , ^ , -.0 ^ ^, ^ ^ , ^ 
L smi L smi 

the triaxiality parameter e being rendered by (!23|) . In fll24p and (11261) . the first two terms 
coincide with those given by the first expression of (2.6) in Kinoshita (1977). They constitute 
^(^jnert) ^ ^j^jj^ ^-^Q third term is /i* . 
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